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There have been considerable interests in the study of low dimensional integrable models 
p|^l6| . One of them is the Hubbard model with 1/r hopping and U on-site energy, introduced 
by Gebhard and Ruckenstein ||^. In this work we report some new results in our recent study 
of the system. 

The long range Hubbard model is exactly solvable in one dimension for arbitrary on- 
site energy [^. Based on small size numerical work and several limiting cases, Gebhard 
and Ruckenstein have introduced an effective Hamiltonian, which was conjectured to be 
equivalent in describing the Hubbard model for arbitrary on-site energy. With the effective 
Hamiltonian, the system has been found to exhibit a metal-insulator transition when the 
bandwidth t is equal to the on-site energy U in the half filling case, and the Gutzwiller state is 
an eigen-state of the effective Hamiltonian in the large U limit. The complete wavefunctions 
and the integrability conditions have remained unsolved for a long time. 

In this primary report, we only concentrate on the study of the wavefunctions and the 
spectrum in the strong interaction limit U = oo. In this extreme limit, one unfortunate 
artifact is that the spins are completely free and decoupled trivially from the charge degrees 
of freedom. We introduce a set of generalized Gutzwiller- Jastrow wavefunctions, and we 
show that they are exact eigen-functions of the Hamiltonian. The lowest energy state in 
this set is also a ground state of the system in the whole Hilbert space. Considering more 
generalized Jastrow wavefunctions, we are able to write the full spectrum. The separation 
of spin and charge in the full excitation spectrum shows that the system is a special example 
of Luttinger liquids in the sense of Haldane. 

The Hamiltonian for the one- dimensional Hubbard model is given by 

H= tijctcja + UY^ni^nn, (1) 

where c}^ and Cjo- are creation and annihilation operators at site i with spin component a. 
We take tij = it(—l)^^^^^/d(i — j) where d{n) = ^sm{mr/L) is the chord distance 0. Here 
we assume periodic boundary condition for the wavefunctions for odd L, or anti-periodic 
boundary condition for even L. 
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In the strong interaction limit U = oo, each site can be occupied by at most one electron. 
In this limit, the Hamiltonian can be written in terms of the Hubbard operator as follows: 

H= i: U.XfXf. (2) 

«7^j;o"=T>-l 

In the following, where-ever in case of the strong interaction, we always implicitly assume 
no double occupancy. Let us denote the number of holes by Q, the number of down-spins 
by M. Following notations used in previous literatures, the state vectors can be represented 
by creating spin and charge excitations from the fully polarized up-spin state \P > [|T0[] , 

1$ >= y: <f (k), ivA) miu h]\p >, (3) 

a,j a j 

where b}^ = cl^iCai is the operator to creat a down-spin at site a, and hj = Cj^ creats a 
hole at site j. The amplitude {Vj}) is symmetric in the down-spin positions, and 

antisymmetric in hole positions. 

To describe uniform motion and magnetization, we consider the following generalized 
Gutzwiller-Jastrow wavef unctions, 

27ri 

^{x, y; Js, Jh) = exp —{Js 'Yxa + JhYl Vi) ^ 

^ a i 

*o = n '^^(^« ~ -Yid^xa - ■ n ^(yi ~ (4) 

a<f3 ai i<j 

where the quantum numbers Jg and Jh govern the momenta of the down-spins and holes, 
respectively. They can be integers or half integers so that we have appropriate periodicities 
(or anti-periodicities) for the wavefunctions under the translations — + L, or i/i — >■ 
Hi + L for odd L (or even L). For the wavefunctions to be exact eigen-functions of the 
Hamiltonian, the momenta of down-spins and holes must take values from some restricted 
regions, which will be specified below. 

These Gutzwiller-Jastrow wavefunctions have been studied extensively in the recent 
study of the integrable models of long range interaction P,[7|JT0HT2| . They were found to be 
the exact eigen-functions of those systems. In the following, we demonstrate that they are 
also the exact solutions of the Hubbard model. To show that they are eigen-functions of the 
Hamiltonian, we have to treat the hopping operator very carefully. 
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The hopping operator can be broken up into two parts, the up-spin hopping operator 
and the down-spin hopping operator. For the up-spin hopping operator = tij'^i'^jh 
its effect is just the hopping of holes alone when it operates on the wavef unctions. But 
the down-spin hopping operator, T| = J^ij^jtijcliCji, will involve the hopping of down-spins 
and holes simultaneously, which needs to be treated using the spin-rotated version of the 
Gutzwiller-Jastrow wavef unctions, developed in the recent work on the supersymmetric t-J 
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For the up-spin hopping, we have 

T^^ix,y;Js,Jh) y.^ (-1)" nj, TT ^.j^) TT F (n) (5) 

where z = exp Fij{n) = cos^nn/L) + sm^nn / L) cot 9ij, Fi^in) = cos^nn/L) + 
sm{7m/ L) cot 9 ia, with % = 7r(?/j — yj)/L,9ia = niyi — Xa)/L. The sum can be calcu- 
lated in a standard way, by expanding the products and classifying the terms by the number 
of particles involved. The terms with more than two particles vanish, yielding the following 
result, 

rt$(x,y; J., J,) ^ _(2^^/^)^^g ^ (27rt/L)^ (6) 
^{x,y;Js,Jh) 

where the hole momenta satisfy the condition \ Jh\ < L/2 — (M -|- Q)/2. 

To deal with the down-spin hopping operators, using the spin rotated version of the 
Gutzwiller-Jastrow wavef unctions, we find 

where Fi^{n) = cos{Txn/L) + sin(7™/L) cot6'j^, Jh = Jh - Js + L/2, 9i^ = Tx{yi - u^)/L, for 
fi = 1,2, ■ ■ ■ , L — M — Q. Here Ui,U2, - ■ ■ ,ul-m~q are the positions of the up-spins on the 
lattice. Finally, we obtain 

^i^^'^'^'/'V'^^^ = -(27rt/L) J,g + i2nt/L)t ^ cot 9,,. (8) 
Here |4| < L/2 - {M + Q)/2, M = L - Q - M. 



After adding the spin-up and spin-down hopping operator effects together, the two par- 
ticle terms vanish due to the fact that the down-spin electrons, the up-spin electrons and the 
holes span the entire lattice. As a result, we see that the Gutzwiller-Jastrow wavefunctions 
are eigenstates of the Hamiltonian, with eigen-energies given by 

E{Js, Jh) = -{'^T^t/L)[2Jh -Js + L/2]Q. (9) 

If we take t to be positive, the lowest energy in this set of Jastrow wavefunctions is 
obtained when Jh and Jh reach their allowed maximum possible values, which occurs at 
Jh = L/2 — (M + Q)/2 and Jg = L — M — Q/2. Thus the dependence of ground state energy 
on the number of holes Q is given by: 

Eo^-{27rt/L)[L/2-Q/2]Q. (10) 

This energy is also the absolute ground state energy of the system in the whole Hilbert 
space, as indicated by our numerical results of 6 site and 8 site lattices. For this ground 
state wavefunction, the exponents of the long distance behaviors of various correlators can be 
obtained from the recent work by Kawakami. In general, the ground state is not unique, and 
this Jastrow wavefunction ground state is only one of them. The ground state degeneracy 
is given by (M + M)!/(M!M!). 

As seen above, the energy consists of decoupled down-spin contribution and the hole 
contribution. In principle, we may find more generalized Jastrow wavefunctions as eigen- 
functions of the Hamiltonian. As will be seen below, in our case, we may even write the 
entire spectrum of the system in the full Hilbert space in terms of a set of Jastrow functions. 

To study other excitations, let us consider more generalized Jastrow wavefunctions in 
the following form 

</. = 0,(X,F)0;,(F)x$o. (11) 

Here the functions (pg and 0^ are polynomials of X = {exp(2niXa/L)}, Y = {exp{2niyi/ L)}. 
They are totally symmetric in their arguments, respectively. The degrees of these polynomi- 
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als must satisfy some specific conditions, so that many particle terms vanish when applying 
the hopping operators. The eigen-energy equation thus reduces to 

Q Q 

- {2T:t/L)[Y,d,{4>s4>h) + Y.Udr + L/2)4>h] = Ect>sct>h. (12) 

i=l 1=1 

where di — YidjdYi. This eigen- value equation can be solved exactly, yielding the spectrum 
given by 

Q Q 

E = -{2Tit/L)[Y,n, + Y,m^]- (13) 

Here, the integers (or half integers ) satisfy the conditions |nj| < L/2 — (M + (5)/2, |m^| < 
L/2 — (M + (5)/2, where rii < rij+i and < m^+i. This result shows that the spectrum 
is invariant when changing the sign of t. 

We may write the spectrum in terms of a set of conjugate quantum numbers defined by 
Ki — Ui -\- rrii -\- {L — Q) /2 + i. The spectrum is given by 

Q 

E = -(27ri/L) {ntQ/L){L + 1), (14) 

1=1 

where Ki takes values from (1, 2, • • • , L). They may be regarded as the momenta of quasi- 
particles "holons" . Our numerical results of 6 site and 8 site lattices indicate that the above 
spectrum spans the full spectrum in the entire Hilbert space. 

In the spectrum formula, each energy level is determined by a charge configuration, such 
as (101010) for L = 6, Q = 3, where the I's represent the charge momenta. We may regard 
the empty values as the momenta for the quasi-particles of spin degrees. Our numerical 
result shows that the degeneracy of each energy level is given by the number of the ways to 
distribute L — Q — M spins s = +| and M spins s — among the empty values. Thus 
we see that the spin degree is decoupled from the charge degree of freedom for the entire 
excitation spectrum. 

For fixed number of electrons Ne — L — Q on the lattice, the free energy consists of two 
parts, F = Fi — TiVg In 2, where the second term comes from the decoupled spin degrees 
of freedom, and Fi is the contribution from the charge degree of freedom, which is that 
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of Q spinless fermions with the relativistic spectrum. In the grand canonical case p7|-[T9 
denoting the chemical potential of the electrons by /i, we find that the free energy per lattice 
site is given by 

F(r, = -/i - — r dq ln[2 + e^^^'^^'\ (15) 
27r J-TT 

Here the number of electrons as a function of the chemical potential is given by N^, = 
— We have also checked that the free energy correctly reproduces the first three 

terms in the high temperature perturbation expansion. 

We have seen that the spin degrees of freedom decouple from the charge degrees of 
freedom for this system in the full excitations in the strong interaction limit. One similar 
system has been studied by one of the authors. He has introduced another integrable model, 
the multi-component Hubbard model of 1/r hopping and U on-site energy, which is the direct 
generalization of the model to the SU{N) case, where a set of generalized SU{N) Gutzwiller- 
Jastrow wavefunctions has been shown to be exact eigen-f unctions of the Hamiltonian in 



the strong interaction limit ]r5|. 



In summary, we have studied the wavefunctions of the one dimensional 1/r Hubbard 
model in the limit of strong interaction U = oo. We have shown that a set of Gutzwiller- 
Jastrow wavefunctions are the eigenstates of the Hamiltonian, and the lowest energy in this 
set is the ground state energy in the whole Hilbert space. The full spectrum can be written 
in terms of more generalized Jastrow wavefunctions. The spin-charge separation occurs in 
the full excitation spectrum, and the system may be regarded as a special Luttinger liquid. 

Finally, we would like to stress that the model is completely integrable for arbitrary 
on-site energy U 0. The most interesting thing is to investigate the wavefunctions and 
the integrability conditions for this 1/r Hubbard model at finite on-site energy U. It would 
be also of great interest to study the ground state properties, such as the spin and charge 
susceptibilities, and various ground state correlation exponents dependence on the on-site 
energy U. For these finite on-site energy studies, further investigations are needed. 
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